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ABSTRACT 
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The  hydroelastic  model  of  the  cochlea  is  used  t.o  analyse  the 
transient  response  to  a  sound  signal.  Using  the  relatively  high 
frequencies  in  the  hearing  range  the  equations  are  first  reduced 
using  viscous  boundary  layer  theory.  The  reduced  problem  is 
solved  by  first  applying  the  Fourier  transform  and  then  using 
the  solution  for  a  pure-tone  signal .  The  response  to  a  tone-pip 
is  then  studied  and  comparisons  are  made,  when  possible,  with 
experiment.  For  example,  the  ’‘ringing'1  observed  in  the 
displacement  of  the  basilar  membrane  is  seen  in  the  model  and 
the  theory  agrees  reasonably  well  with  the  neural  latency 


measurements . 
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INTRODUCE IJN 

In  the  study  o-f  the  transmission  and  resolution  of  the  sound 
signal  in  the  cochlea  the  behavior  associated  with  the  longtime 
response  to  a  pure  tone  has  received,  by  far,  the  most  attention. 
This  is  understandabl e  as  the  work  on  this  aspect  of  the  hearing 
process  has  been  invaluable  in  the  elucidation  of  the  basic 
mechanisms  underlying  frequency  discrimination.  However,  to  be 
able  to  address  questions  connected  with  speech  discrimination  and 
noise  induced  hearing  loss  it  is  necessary  to  consider  the 
transient  behavior.  Unfortunately,  the  effort  that  has  been 
invested  in  this  aspect  of  the  hearing  process  has  been  limited. 

In  terms  of  direct  experimental  observation  of  the  transient 
response  of  the  basilar  membrane,  Robles,  et  al  (1976)  is  one  of 
the  only  studies  presently  available.  In  fact,  most  of  the 
experimental  works  on  the  transient  problem  give  only  indirect 
indications  of  the  motion  of  the  basilar  membrane.  This  includes 
psychoacoust i cal  tests,  such  as  Stapells,  et  al  (1982),  and  neural 
measurements ,  such  as  Anderson,  et  al  (1970)  and  Pfeiffer  and  Kim 
(1972).  The  modeling  efforts  are  also  limited  and  the  few  that 
have  appeared  in  the  literature  have  concentrated  almost 
exclusively  on  the  low  frequency  response.  In  fact,  there  has 
never  been  a  complete  analysis  of  the  fully  three  dimensional 
transient  motion  in  the  cochlea. 

It  is  for  these  reasons  that  the  objective  of  the  present 
study  is  to  develop  a  description  of  the  transient  behavior  in  the 
cochlea  that  covers  the  entire  hearing  range.  The  model  that  is 


used  for  the  cochlea  is  the  same  as  the  one  that  has  been  used  to 


describe  the  response  to  a  pure  tone  (Holmes  and  Cole,  1983,1984; 
Holmes,  1986).  The  formulation  of  the  model  is  relatively  simple 
and  it  assumes  the  fluid  filling  the  chambers  is  viscous  and 
incompressible  and  the  basilar  membrane  is  modeled  as  a  fibrous 
elastic  plate.  Except  for  the  fluid's  density  and  viscosity,  and 
the  transverse  elastic  modulus  and  mass  density  for  the  basilar 
membrane,  all  of  the  parameters  in  the  model  are  geometric.  The 
values  for  these  four  parameters  are  known  from  independent  tests 
so  the  model  of  the  mammalian  cochlea  used  here  is  characterized 
entirely  by  the  geometric  structure  of  the  ear.  Because  of  this 
the  model  can  be  applied  easily  to  a  wide  variety  of  mammalian 
cochlae  (Holmes,  1986). 

In  the  next  section  the  equations  governing  the  motion  in  the 
cochlea  are  presented  and  nondi mensi onal i zed .  The  reason  for 
scaling  the  problem  is  to  simplify  the  analysis  used  to  solve  it. 
This  is  done  by  first  using  viscous  boundary  layer  theory  to 
reduce  the  problem.  Once  this  is  done  a  slender  body 
approx i mat i on ,  which  is  based  on  the  long  narrow  geometry  of  the 
cochlea,  is  used  to  obtain  the  solution  of  the  problem.  This  is 
the  same  procedure  that  was  used  to  study  the  response  to  a  pure 
tone  although  the  details  of  the  derivation  are  somewhat  different 
due  to  the  time  dependence  of  the  problem.  In  any  case,  once  the 
solution  is  derived  it  is  then  a  simple  matter  to  analyze  to 
response  of  the  basilar  membrane  to  various  signals. 

The  reason  that  viscous  boundary  1 ayer  theory  can  be  used  is 
because  of  the  relative  thinness  of  the  layer  as  compared  to  the 
cross-sectional  dimension  of  the  fluid  chamber.  This  should  not  be 


interpreted  to  mean  that  the  fluid  viscosity  is  unimportant.  In 
the  model  considered  here  the  fluid  viscosity  is  the  only 
dissipative  mechanism,  and  with  it  the  theory  gives  a  very 
accurate  description  of  the  frequency  map,  the  moderate  intensity 
amplitude  and  phase  tuning  curves,  and  the  static  stiffness 
measurements  (Holmes,  1986).  This  includes  human  as  well  as  other 
mammals  (e.g.,  the  cat,  guinea  pig,  and  chinchilla).  At  the  same 
time,  with  the  exceptions  noted  earlier,  the  parameters  that 
define  the  model  are  all  geometric  and,  as  shown  in  Holmes  (1985), 
they  are  determined  solely  from  the  measured  values  for  each  of 
these  animals.  This  raises  the  question  as  to  the  importance  of 
other  passible  dissipation  mechanisms,  such  as  membrane  damping. 
Whatever  contribution  they  make  it  is  in  addition  to  fluid 
viscosity.  However,  one  of  the  difficulties  with  attempting  to 
account  for  other  damping  mechanisms  is  that  there  are  no  direct 
experimental  measurements  on  which  to  base  the  modeling 
assumptions.  For  example,  there  has  never  been  a  study  of  mode 
separation  and  transient  amplitude  decay  for  the  basilar  membrane. 
This  makes  the  modeling  effort  almost  guesswork,  although  one  way 
around  this  is  to  assume  that  the  damping  mechanism  is  similar  to 
that  found  in  other  fibrous  biological  membranes,  such  as  the 
tympanic  membrane  (Rabbitt  and  Holmes,  1985).  Except  for  Allen 
(1980)  this  has  never  been  done  and  the  usual  assumption  is  that 
the  damping  is  simply  a  dashpot  mechanism.  Such  an  assumption  is 
difficult  to  justify  as  it  has  little  basis  in  the  ul trastructure 
of  the  basilar  membrane  or  the  organ  of  Corti . 


EQUATIONS  OF  MOTION 

The  assumptions  on  the  geometrical  and  structural  components 
of  the  cochlea  are  similar  to  those  made  in  Holmes  and  Cole 
(1984).  Consequently,  the  presentation  to  follow  is  relatively 
brief  but  a  more  complete  discussion  of  the  assumptions  made  here 
can  be  found  in  the  earlier  paper.  The  cochlea  is  taken  to  be  a 
tapered  tube  containing  two  fluid  filled  chambers  (Fig.  1).  The 
partition  separating  these  chambers  consists  of  a  rigid  section 
(the  bony  shelf),  a  flexible  portion  (the  basilar  membrane)  and  an 
aperture  at  the  apical  end  (the  hel i cotrema) .  The  cochlear  wall 
is  rigid  except  at  the  basal  end  where  there  are  two  openings  (the 
oval  and  round  windows)  that  are  covered  by  membranes.  The 
footplate  of  the  stapes  is  attached  to  the  cochlea  at  the  oval 
window.  The  variation  of  the  geometry  is  fairly  arbitrary  except 
that  it  is  assumed  that  the  partition  and  basal  end  are  planar 
surfaces  and  the  cochlear  wall  is  symmetric  through  the  partition. 

The  fluid  filling  the  two  chambers  is  similar  in  terms  of  its 
mechanical  behavior  to  water  at  body  temperature.  Accordingly,  it 
is  assumed  to  be  an  incompressible  Newtonian  fluid,  and  so,  the 
equations  of  motion  are 

K*  -  \  vp*  .  <>•’ 

and 

V-  V*  =  0  ,  (lb) 

where  v*(x*,  y*,  z*,  t*)  is  the  velocity,  p#(x*,  y*,  z*,  t*> 
is  the  pressure,  r  is  the  kinematic  viscosity  and  p  is  the  mass 
density.  The  asterisk  appearing  in  (1)  is  used  to  distinguish  the 
variables  from  their  nondi mensi onal  analogs  introduced  below. 
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OVAL  WINDOW 


The  geometry  and  notation  used  for  the  hydroelastic  model  of  the  cochlea 


Also,  the  notation  used  here  is  the  same  as  in  Holmes  and  Cole 


(1984)  except  for  the  displacement  of  the  stapes  which  is  given  in 
(4)  below. 

The  basilar  membrane  (BM)  is  modeled  as  an  inhomogeneous, 
strongly  orthotropic,  elastic  plate.  In  this  case  the  equation  of 
motion  is  (Lekhni tski i ,  1968) 


D*«Jy*n*  +  U*3^**)*  =  ~  [£  p*  (x*,y#,0,t*)  JJ 


where  t)*(x  ,  y  ,  t* )  is  the  displacement  in  the  z*-direction ,  D; 

is  the  transverse  bending  rigidity,  and  n*  is  the  mass  per  unit 
area  of  the  basilar  membrane.  The  longitudinal  bending  and 
twisting  terms  have  not  been  included  in  (lc)  as  they  have  been 
found  to  be  negligible  in  the  basilar  membrane  (Voldrich,  1978; 
Holmes  and  Cole,  1984).  The  right-hand  side  of  (lc)  represents 
the  fluid  loading  on  the  BM  and  is  given  as 


=  p* (x* , y* ,0+ , t*)  -  p*(x*,y*,0  ,t*> 


By  assuming  the  BM  is  inhomogeneous  it  is  meant  in  (lc)  that  its 
material  coefficients  depend  on  the  longitudinal  variable  x*. 
Also  note  that  the  displacement  of  the  basilar  membrane  is 
positive  if  it  is  in  the  positive  z *-di recti  on. 

To  complete  the  description  of  the  problem  the  boundary  and 
initial  conditions  need  to  be  specified.  Because  the  fluid 
velocity  must  be  continuous  at  the  boundary  we  have  that 


3*  = 


on  rigid  wall 


(0,0,dt*n*)  on  BM 


"V7V  .VfcV.V  V  V  A  V  \ 


Similarly,  at  the  aval  and  round  windows 

3*  =  o,  0)  ,  <4> 

where  n|(t#)  represents  the  displacement  of  the  stapes  footplate. 
In  (4),  Us  is  positive  if  it  is  in  the  positive  x*-direction.  For 
the  boundary  conditions  for  the  plate  it  is  assumed  that  it  is 
simply  supported.  It  should  be  pointed  out,  however,  that  the 
methods  used  to  solve  the  above  problem  do  not  depend  on  the 
particular  boundary  conditions  for  the  basilar  membrane  and  one 
could  just  as  well  assume  it  is  clamped  or  cl amped/si mpl y 
supported  (Holmes,  1986).  As  for  initial  conditions,  it  is 
assumed  that  the  system  starts  from  rest  from  its  equilibrium 
posi t i on . 

The  transient  problem  for  the  cochlea  is  now  complete  and  it 
involves  solving  the  fluid  equations  (la,b)  coupled  with  the 
elastic  boundary  (lc).  Given  the  three  dimensional  geometry  and 
the  relatively  short  waves  that  propagate  along  the  basilar 
membrane  even  a  numerical  solution  of  this  problem  would  be 
difficult  as  well  as  impractical  at  this  time.  However,  the 
cochlea  has  a  long  slender  geometry  and  the  audible  frequency 
range  is  relatively  high.  Both  of  these  can  be  used  to  reduce  the 
problem  but  to  do  so  it  is  necessary  to  scale  it.  As  in  Holmes 
and  Cole  (1984)  the  spatial  coordinates  are  nondi mensional ized  as 
foil ows 

x*  =  Lx,  y*  =  By,  and  z*  =  Bz , 
where  B,L  are  the  width  and  length  of  the  basilar  membrane, 
respectively.  Also,  setting  t*  =  tct. 
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(5) 


where  n50  is  the  amplitude  of  the  stapes  and  uc  is  a 
character i stic  value  -for  u*  (e.g.,  it's  value  at  x=l/2).  The 
characteristic  time  scale  tc  used  here  is 


B\ic 

D*c 


(6) 


where  D2c  is  a  characteristic  value  of  D*.  Thus,  tc  is  a  measure 
of  the  fundamental  frequency  of  an  elastic  beam  in  a  transverse 
cross-section  of  the  basilar  membrane.  For  example,  for  the  human 
ear  tc  *v  1.3  msec.  In  Holmes  and  Cole  (1984)  the  characteristic 
time  scale  was  the  driving  frequency  of  the  stapes.  Although  this 
differs  from  (6)  the  results  of  the  analysis  to  follow  actually 
contains  this  earlier  work. 

In  dimensionless  coordinates  the  problem  describing  the 
motion  in  the  cochlea  is: 
i )  fluid 


“  S2VfjO  =  -  Vp  ,  (7a) 

£  Z®xvl  +  ®yV2  +  ®zv3  =0  ,  (7b) 


ii  > 


basi lar 

DjdylJ  + 


membrane 
,  * 

U  3  fcT1 


,y»o,t) 


11 


(8) 
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iii)  boundary  conditions 


-» 

v  = 


<0,0,aatn> 

<att>s,o,o) 


on  rigid  wall  (9a) 
on  BM  (9b) 
on  stapes  footplate  .  (9c) 


The  initial  conditions  to  be  used  are  that  the  motion  starts 


rest  at  t=0.  In  (7a)  ,  V*  =  £  23*  +  + 

D2  =  D*/D2c  anc*  U  =  U*/Uc-  Also, 


In  (B), 


6  2  = 


-c 

>2  ’ 


L  * 


and  a  = 


The  last  parameter  is  simply  the  ratio  of  the  mass  of  the  fluid 
above  the  basilar  membrane  to  the  characteristic  mass  of  the 
plate.  As  for  £  ,  it  is  the  ratio  of  the  width  to  length  of  the  BM 
and  it  is  relatively  small  (e.g.,  for  humans  £  lO- 2 ) .  This  is 
used  later  to  derive  a  slender  body  approximation  to  the  solution 
of  the  problem.  The  parameter  8  is  the  inverse  of  the  Reynolds 
number,  where  the  latter  is  based  on  the  cross-sectional  width  of 
the  chamber.  This  too  is  relatively  small  and,  accordingly, 
viscous  boundary  layer  theory  will  be  used  to  reduce  the  above 
problem. 

The  reduction  of  the  transient  problem  consists  of  two 
steps.  First,  viscous  boundary  layer  theory  is  used  to  reduce  the 
fluid  problem.  From  this,  rather  than  solving  (7),  we  only  need 
to  solve  Laplace's  equation  for  the  pressure,  where  the  viscosity 
contributes  through  a  modified  boundary  condition.  This  reduced 
problem  is  then  solved  using  a  slender  body  approximation  based  on 
the  fact  that 


c 


BOUNDARY  LAYER  APPROXIMATION 


The  frequencies  in  the  hearing  range  are  fairly  high,  for 
example,  they  range  from  about  20  to  20,000  Hz  for  the  human  ear. 
This  serves  as  the  basis  for  the  first  reduction  of  the  equations 
of  motion  and  involves  the  application  of  viscous  boundary  layer 
theory.  In  the  approximations  to  follow  it  is  assumed  that 
8  <<  1.  Also,  in  the  discussion  the  other  parameters  in  the 
problem  are  assumed  fixed.  The  analysis  splits  naturally  into 
three  parts,  the  first  two  involve  deriving  the  approximations  for 
the  inviscid  and  boundary  layer  regions  and  the  third  step  deals 
with  their  matching.  The  details  are  somewhat  involved  but  they 
result  in  a  considerable  reduction  of  the  fluid  equations.  For 
those  who  may  wish  to  skip  the  derivation,  the  conclusions  of  the 
analysis  are  given  in  the  last  paragraph  of  this  section, 
i)  Inviscid  Region:  In  the  inviscid  core  the  appropriate 
expansions  of  the  fluid  velocity  and  pressure  are 

v  ~  v0(x,y,z,t>  +  8vt  +  ...  ,  (10a) 

and 

p~p0(x,y,z,t)+6p1+...  .  (10b) 

Substituting  these  into  (7)  and  equating  like  powers  of  8 
leads  to  the  fallowing  problem  far  the  nth  term  in  the  expansion 

0  ( 6 n >  V*pn  =  0 

H'Jn  =  “Vpn  .  (ID 

As  expected  it  is,  in  general,  not  possible  to  satisfy  all  of  the 
boundary  conditions  with  this  approximation  (in  particular,  the  no 
slip  condition  in  (9)  can  not  be  satisfied).  This  necessitates 


the  introduction  of  a  boundary  layer. 

ii)  Boundary  Layer:  To  obtain  the  boundary  layer  approximation 
only  the  layer  above  the  partition  is  considered.  The  analysis 
for  the  rest  of  the  cochlear  wall  is  similar  and  involves  the 
introduction  of  orthogonal  coordinates  orientated  with  respect  to 
the  normal  direction  (Holmes  and  Cole,  1984).  For  the  region 
above  the  partition  the  boundary  layer  coordinate  is 

z  =  /s  .  (12) 

The  Navi er-Stokes  equations  (7)  in  this  case  take  the  form 

[«  ”  ”  8  2  V  * ,  y)  v  =  -(«xi«y.f  «~]p  ♦  (13a> 

d^v,  +  S(C2®xv1  *  9yV,)  =  0  (13b) 

where 

Vx , y  «  « 2«x  +  - 

The  appropriate  expansions  in  this  layer  are 

v  'v  0Q(x,y,z,t)  +  S\Jt  +  ...  ,  (14a) 

and 

p  ~  P0(x,y,z,t)  +  6Pt  +  _  .  (14b) 

Proceeding  as  before,  by  substituting  these  into  (13)  and  equating 
like  powers  of  8  the  following  problems  are  obtained: 

0(1)  3~P0  =  O 

«-VJ0  =  0 

(*t  ~  ,2’)V'S0  =  “  *sPo  ’  <s  =*’V> 

Solving  these  equations  and  using  boundary  condition  (9b)  one 


finds  that 


=  Po <*  »y»t> 


V 


Vj o  ~  “ o  » 


®sP0dT  ,  <s=x,y) 


To  obtain  a  uniformly  valid  two  term  approximation  it  is 
necessary  to  solve  for  at  least  some  of  the  second  terms  in  the 
expansion. 


□  <6  ) 


i  =  -  “3. Ho 


v*-y 


Po  dr 


Using  the  Laplace  transform  one  finds  that  the  solutions  of 
the  above  problems  are 

Pt  =  -oz3*ti0  +  b(x,y,t) 


Vji  =  +1  jz  -  2  v  t-T  — i  erf  c 


[stM]  v“- 


^odT  * 


where  b  is  a  constant  of  integration. 

iii)  Matching:  It  remains  to  match  the  expansions  from  the  two 
regions  and  to  do  so  an  intermediate  variable  is  used,  which  is 
defined  as 


z  V  S 


With  this,  using  Taylor's  theorem  the  inviscid  approximation  for 
the  pressure  (10b)  can  be  written  as 


P"  Po  lz*o+  '/^-,Z*zPo  lz*o+  6  (Pl  +  2  2*dzPo>  |2=0  +  •  *  •  »  <16> 


and  the  boundary  layer  approximation  (14b)  is 

p  **•  P0  -  oz/T‘«*n0  +  Sb  +  ...  .  (17) 

Matching  (16)  and  (17)  as  well  as  matching  the  velocities  one  finds 


P°  "  Polz=o 


dzPo  lz=< 


=  » 


*zPllz=o  ST*  J*0  Vx’ 


y  Po  dT 


Combining  these 


9’P‘U  '  -  J  yf=? 


Po (x,y,0,t)  dr  .  (18) 


Based  on  the  above  analysis,  a  uniformly  valid  two  term 
approximation  (in  8)  of  the  fluid  pressure  is  obtained  by 


solving 


Vf  P  =  °  » 


where,  for  O  <  x  <  1, 


«  p  -  ~=7 

n  V  tt ' 


f  =L 

J  y  t-n 

*  n 


dnpdr  = 


on  rigid  wall 


(20a> 


on  BM 


In  (20a) ,  n  is  the  unit  outward  normal  to  the  cochlear  wall.  A 
similar  condition  applies  at  the  basal  end  (x=0)  and  is  given  as 


*xP  +  v^  I  3xpdT  =  "  V*  •  (20b> 

*  0 

Therefore,  the  original  equations  for  the  fluid  (7)  along  with 
boundary  conditions  (9)  have  been  replaced  with  (19)  and  (20).  In 
other  words,  it  remains  to  only  solve  one  equation  for  the  fluid 
and  the  viscosity  now  only  contributes  in  the  boundary  condition, 
rather  than  in  the  equation  of  motion.  It  is  also  of  interest  to 
note  that  this  approximation  reduces  to  the  one  derived  earlier  by 
Holmes  and  Cole  (1984)  for  the  longtime  periodic  response  to  a 
pure— tone  forcing  as  t  -»  «.  Moreover,  it  reduces  to  the  transient 
low  frequency  approximation  used  in  Holmes  (1981,  1982)  in  the 
case  of  a  low  frequency  signal. 


SOLUTION  OF  REDUCED  PROBLEM 


It  remains  to  solve  the  reduced  fluid  equations  derived  in 
the  last  section  in  conjunction  with  the  plate  equation  and 
associated  boundary  conditions.  There  are  a  number  of  ways  to  do 
this  and  the  one  to  be  used  here  involves  the  Fourier  transform  in 
t.  For  a  variable  such  as  n(x,y,t)  its  transform,  which  is 
represented  with  a  capital  N(x,y,u)  as  well  as  by  the  symbol  dMn), 
is  defined  as 


To  be  able  to  transform  the  equations  of  motion  it  is 
necessary  to  assume  that  the  Fourier  transform  of  the  stapes 
displacement  is  well  defined.  With  this  the  transf ormation  of  the 
problem  is  fairly  straightf orward,  except  perhaps  for  the 
convolution  integral  in  boundary  conditions  (20).  The  latter  can 
be  evaluated  by  interchanging  the  orders  of  integration  as  follows 


o 


(21) 


Vt  P  =  O  , 

and  equation  (8)  for  the  BM  is 

<D2«y  -  u2u ) N  =  — 2P (x,y,0,u)  , 

where,  from  (20a) , 


(22) 


®nP  " 


y/  i(J 


8  *Ap  = 


i  “n' 


0  on  rigid  wall 

.  -«uJN  on  BM  , 


and,  from  (20b) , 

c  8  t  . 

dxP  +  1  ~  =  u*Ns  on  stapes  footplate 

-/  iu ' 

This  problem  is  essentially  the  same  as  the  one  studied  in  Holmes 
and  Cole  (1984).  Consequently,  it  is  possible  to  simply  write 
down  a  first  term  approximation  of  the  solution.  This 
approximation,  which  is  based  on  the  assumption  that  E  <<  1,  is 
derived  using  a  multiple  scale  argument.  The  result  is  that  the 
solution  far  the  (transformed)  BM  displacement  has  the  form 


N  **•  A  ( x  ,  «j)  e  ie<X,U>T|0(x,y)  .  (23) 

The  functions  in  this  expression,  such  as  A(x,u),  will  be  given 
later  as  they  are  not  needed  to  complete  the  derivation  of  the 
solution  of  the  problem.  Inverting  (23)  using  the  convolution 
theorem  leads  to 


'  . — K0  (x  ,y,t  — t>  n, 

y  2w  j 


(  t  )  dr 


(24a) 


where 

K0  <* ,y»t) 


J  1 


-  A  <  x 

a 


.  -i  6  <x  ,  (ii)  | 
,u)e 


n  o <x ,y ) 


(24b) 


Note  that  it  has  been  assumed  in  (24)  that  n  =  0  for  t  ^  O  and 
the  function  ti0  is  independent  of  u. 


fi 


i 


»y 


It  remains  to  invert  the  inverse  Fourier  transform  in  (24b>. 
However,  because  of  the  relatively  complicated  nature  of  the 
functions  there  is  no  apparent  closed  farm  inversion. 
Nevertheless,  there  are  some  simplifications  that  can  be  made. 

For  example,  since  A  is  an  odd  function  of  u  and 

0(x,-u)  =  — 0 (x , w)  +  *  then  it  is  only  necessary  to  integrate  over 
0  <  u  <  ®  in  (24b) . 


k 

i 


* 
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SOLUTION  OF  TRANSIENT  PROBLEM 


In  the  last  two  sections  the  original  equations  of  motion 
were  reduced  and  a  first  term  approximation  to  the  transient 
problem  was  derived.  The  result  is  that  the  displacement  of  the 
basilar  membrane  is  given  as 

CD 

tt*(x,y,t)  J*  K(x,y,t-T)  its<T)dT  ,  (25a) 

o 

where 

K(x,y,t)  =  *  ^  cos  |ut  -  0(x,u)jdu  *  n0(x,y).  (25b) 

o 

The  amplitude  A(x,u)  and  the  phase  6(x,u)  are  given  in  (A18)  and 
(A19)  ,  respectively.  The  function  ti  0 ,  which  is  the  displacement 
of  the  wave  in  the  transverse  cross-sect i on  of  the  BM,  is  given  in 
(A9)  .  The  time  history  of  the  displacement  of  the  BM,  as  given  in 
(25),  is  in  the  form  of  a  convolution  integral.  Because  of  the 
complexity  of  the  functions  in  (25b>  it  appears  to  be  necessary  to 
evaluate  (25)  numerically.  However,  the  kernel  defined  in  (25b) 
needs  to  be  calculated  only  once.  To  do  so  for  a  given  x  and  o 
one  must  find  the  solution  of  a  transcendental  equation  (AlO) , 
called  the  dispersion  relation,  for  the  local  inviscid  wave  number 
k0.  Once  k„(x,u)  is  known  then  it  is  a  relatively  simple  matter 
to  calculate  A  and  0  from  the  formulas  given  in  the  Appendix. 

In  the  calculations  to  follow  for  the  human  ear  the  kernel  is 
evaluated  at  100  equally  spaced  points  along  the  x-axis,  at  50 
equally  spaced  frequencies  between  0  and  500  Hz,  and  at  50 
1 ogar i thmi cal  1 y  spaced  frequencies  between  500  and  15,000  Hz.  The 
methods  used  to  evaluate  the  kernel  are  the  same  as  used  in  Holmes 
and  Cole  (1984)  and  Holmes  (1986).  The  reason  the  frequency 


spacing  is  divided  in  this  way  is  because  of  the  effects  of  the 


low  frequency  modes  (Holmes,  1985).  In  particular,  the  amplitude 
and  phase  functions  are  not  monotonic  for  low  frequencies  so  it  i 
necessary  to  increase  the  number  of  points  in  this  region  so 
interpolated  integration  can  be  used  in  (25b) .  By  doing  this,  it 
takes  about  85  sec  of  cpu  time  on  an  IBM  3081D  to  calculate  the 
amplitude  and  phase  functions  in  (25b).'*' 


+The  FORTRAN  programs,  and  documentation,  used  to  calculate  the 
solution  of  the  transient  problem  are  available  upon  request  from 


the  author. 


APPLICATION  TO  THE  HUMAN  COCHLEA 


The  transient  response  of  the  basilar  membrane  is  now 
illustrated  by  evaluating  (25)  tor  the  human  ear.  To  do  so 
the  values  used  -for  the  -fluid  kinematic  viscosity  and  density 
are  r=0.008  cm2/sec  and  />~1.0  gm/cm3.  For  the  basilar 
membrane  L=3.5  cm,  the  transverse  bending  rigidity  is 

E, (h*) J 


d; 


12  ( l-<r*  ) 


the  mass  density  is 

u*  =  Uo^h*  ’ 

and  the  thickness  is 

h*  =  h£(l— hsx> 

where  dyn/cm2,  h*=0.74xl0~3  cm,  hs=0 . 38 ,  or=  1/2,  and 

u0  =  17.  The  longitudinal  boundaries  o-f  the  BM  are  y*=±BG(x), 
where  B=0 . 05  cm  and 

G(x)  =  —  (5x  +  1)  ,  for  O  <  x  <  1  . 

Also,  the  transverse  cross-section  of  the  cochlea  chamber  is 
rectangular  with  a  constant  area  of  0.01  cm2  and  the  area  of 
the  stapes  footplate  is  0.03  cm2.  These  values  are  the  same 
as  used  when  studying  the  response  to  a  pure  tone  (Holmes, 
1986)  and  are  representat i ve  of  the  measured  values  for  the 
human  ear . 

The  forcing  function  to  be  considered  is  a  tone-pip  in 
which  the  displacement  of  the  stapes  has  the  form 


where  t*  =  2nir/cj*  and  n  is  an  integer.  Inserting  this  into 
(25)  and  integrating  one  finds  that 

a> 

H  (■)  ^ 

tl*(x,0,t)  ~  - — =-  A  (x  ,  u)  cos  fut-6  (x  ,  u)  -  ^<jt0lsinf^  ut0]du 

*  u(u0~u  )  L  ^  J  J 

0  (27) 

where  t0  =  t*/tc  and  u0  =  tcw*.  The  numerical  evaluation  of 

this  integral  in  the  case  of  when  f=1400  Hz  and  n=2  is  shown 

in  Fig.  2  at  four  successive  times.  It  is  seen  in  this 

figure  that  the  initial  response  of  the  BM  is  in  the  form  of 

a  traveling  wave  that  propagates  towards  the  apical  end.  The 

wave  takes  approximately  9  msec  to  transverse  the  BM.  Also 

note  that  it  is  the  low  frequency  components  of  the 

disturbance  that  reach  the  apical  end  and  they  travel  faster, 

and  farther,  than  the  higher  frequency  waves.  This  is  not 

unexpected  given  the  dependence  of  the  phase  velocity  on 

frequency  (Holmes,  19B6) . 

The  time  histories  of  the  displacement  of  the  BM  at 
three  spatial  locations  are  shown  in  Fig.  3  using  the  same 
signal  as  in  Fig.  2.  In  each  case  there  is  a  time  delay 
before  the  BM  begins  to  move.  To  measure  this,  the  latency 
is  defined  as  the  time  until  the  displacement  reaches  its 
first  significant  positive  peak.  This  choice  is  made  because 
it  approximates  the  value  measured  neurally,  although  the 
latter  also  includes  other  sources  of  delay.  In  any  case,  as 
seen  in  Fig.  3  the  time  it  takes  the  wave  to  reach  x=0.25 
(approximately  0.B5  msec)  is  less  than  it  takes  it  to 
propagate  from  x=0.5  to  x=0.75  (which  takes  about  2.9  msec). 


This  is  due  to  the  monotonic  decrease  of  the  phase  velocity 
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with  longitudinal  position  (Holmes,  1986).  Something  else 
that  is  seen  in  this  -figure  is  that  the  BM  oscillates  much 
more  at  the  point  x=0.5  than  at  the  other  two  locations.  The 
reason  for  this  is  because  the  characteri stic  frequency  for 
x=0.5  is  1400  Hz.  If  the  driving  frequency  is  4570  H2 ,  which 
is  the  characteristic  frequency  for  x=0.25  then,  as  seen  in 
Fig.  4,  the  oscillations  occur  at  x=0.25  rather  than  at 
x=0. 5. 

A  contour  plot  for  the  displacement  of  the  basilar 
membrane  is  shown  in  Fig. 5  for  the  same  signal  that  was  used 
for  Figs.  2  and  3.  Only  positive  displacements  are  shown  in 
this  figure  because  of  their  importance  in  the  depolarization 
of  the  hair  cell  (Bell  and  Holmes,  1986).  In  any  case,  the 
wave-like  nature  of  the  response  is  clearly  seen  along  with 
the  ringing  that  occurs  due  to  the  effects  of  the 
character i sti c  frequency.  The  "chain  of  islands,"  which  run 
from  the  lower  left  to  upper  right,  correspond  to  individual 
wave  packets  propagating  down  the  basilar  membrane.  Except 
for  the  lower  two,  these  chains  are  parallel.  Also,  each 
succeeding  chain  is  smaller  than  the  previous  one  but  x=0.5 
remains  at,  or  near,  the  center  of  each.  The  contributions 
of  the  low  frequencies  appear  primarily  in  the  two  lower 
packets  and  even  then  they  are  confined  principally  to  the 
first.  This  is  due  to  the  start  up  of  the  signal  which 
contains  frequencies  throughout  the  spectrum. 

In  the  case  of  a  pure  tone  signal  it  is  possible  to 
determine  the  length  of  time  it  takes  for  the  longtime 
periodic  solution  to  appear.  With  a  forcing  such  as  given  in 


X-AXIS 


Fig.  5.  The  level  contours  for  the  positive  displacement  of 
the  basilar  membrane,  normalized  by  the  amplitude  of  the 
stapes,  in  response  to  a  tone-pip  with  frequency  1400  Hz.  T 
curves  represent  constant  normalized  displacements  of  10,  20 


Fig.  6.  Transient  response  of  the  basilar  membrane  for  a 
pure-tone  signal  with  frequency  1400  Hz,  at  time  a)  0.5,  b> 
2.0,  c)  8.0,  and  d)  16  msec.  Also  shown  is  the  envelope  of 
the  corresponding  longtime  solution. 
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(26)  this  is  the  time  it  takes  the  transients  associated  with 
the  start-up  of  the  signal  to  propagate  down  the  basilar 
membrane  and  decay.  Based  on  the  discussion  in  the  earlier 
paragraphs  this  is  approximately  5  msec  for  the  higher 
frequency  components  and  about  10  msec  for  the  lower 
frequencies.  The  solution  as  calculated  from  (27)  and  the 
envelope  of  the  longtime  solution  are  shown  in  Fig.  6  at  four 
successive  times  in  the  case  of  when  f=1400  Hz  and  n=50.  As 
expected,  at  the  earliest  time  there  are  still  appreciable 
transient  effects,  the  higher  frequency  transients  have 
disappeared  at  the  three  later  times  and  the  lower 
frequencies  have  effectively  decayed  at  the  last  two. 

Before  leaving  this  example  the  method  used  to  calculate 
the  solution  should  be  explained.  With  the  earlier 
evaluation  of  the  amplitude  and  phase  functions  in  (25b)  the 
integral  in  (25a)  is  evaluated  using  Simpson's  rule  with  400 
equally  spaced  points  between  O  and  15,000  Hz.  This  is  done 
using  the  values  for  the  amplitude  and  phase  functions 
calculated  earlier,  which  were  evaluated  at  a  smaller  number 
of  frequency  points,  by  using  linear  i nterpol ati on .  In  doing 
this  it  takes  approximately  10  sec  of  cpu  time  to  calculate 
all  of  the  curves  shown  in  Fig.  3. 


DISCUSSION 


The  method  that  has  been  used  here  to  solve  the  transient 
problem  is  based  on  the  spectral  decomposition  of  the  sound 
signal  using  the  Fourier  transform.  By  doing  this  the 
transformed  problem  is  essentially  the  same  as  the  one  for 
the  longtime  response  to  a  pure— tone,  so,  the  solution  is 
relatively  easy  to  obtain.  The  factor  that  complicated  the 
analysis  was  the  viscous  boundary  layer  as  it  was  necessary 
to  derive  an  approximation  that  described  the  development  of 
the  layer  in  time.  Once  this  has  been  done  the  solution  to 
the  transient  problem  is  in  the  form  of  a  convolution 
integral  in  which  the  kernel  is  determined  completely  from 
the  response  for  a  pure-tone. 

Although  the  analysis  used  to  obtain  the  solution  is 
straightforward  it  is  necessary  to  evaluate  the  integral 
numerically.  By  comparison,  for  the  low  frequency  theory  the 
solution  is  simple  enough  that  it  is  relatively  easy  to 
derive  the  characteri sties  of  the  response  of  the  basilar 
membrane  (Holmes,  1981).  Nevertheless,  the  calculations  used 
to  evaluate  (25)  are  not  particularly  difficult  and  lead  to 
some  interesting  conclusions.  For  example,  for  the  human 
cochlea  the  time  scale  associated  with  the  development  of  a 
pure  tone  is  about  10  msec.  The  decay  of  the  displacement  of 
the  basilar  membrane  is  also  strongly  dependent  on  the 
spectrum  of  the  signal.  Even  in  the  case  of  a  pure  tone-pip 
lasting  two  cycles  the  basilar  membrane  at  the  characteri st i c 
place  vibrates  through  10  or  15  cycles.  Although  there  are 
no  direct  experimental  measurements  to  compare  with. 


responses  similar  to  the  ringing  seen  in  the  model  have  been 
observed  both  in  the  displacement  of  the  BM  as  well  as  in  the 
neural  response  (Robles,  et  al  ,  1976;  Pfeiffer  and  Kim, 

1972).  Robles,  et  al  also  observed  that  the  latency,  or  time 
delay,  of  the  response  to  reach  x=0.  1  in  the  squirrel  monkey 
is  about  0.35  msec.  The  model  predicts  that  at  the  same 
paint  for  the  human  the  delay  is  only  about  0.06  msec.  It  is 
not  clear  why  there  is  so  much  difference  between  the  two. 

The  theory  also  predicts  that  the  total  time  it  takes  the 
response  to  reach  the  apical  end  is  about  7.75  msec.  From 
neural  measurements  the  travel  time  appears  to  be  between  7 
and  S  msec  for  the  squirrel  monkey  (Anderson,  et  al ,  1971), 

for  the  chinchilla  (Siegel,  et  al ,  19B2) ,  and  for  the  cat 

(Kim  and  Molnar,  1979).  So,  for  this  the  model  is  in 
reasonable  agreement  with  experiment.  It  remains  though  to 
either  extend  the  model  to  account  for  the  hair — cell  system, 
or  for  direct  experimental  observations  of  the  transient 
motion  of  the  BM,  before  a  more  complete  comparison  between 
theory  and  experiment  can  be  made. 
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APPENDIX.  EVALUATION  OF  THE  KERNEL  FUNCTION 

Once  the  transient  problem  is  reduced  using  viscous  boundary 
layer  theory,  and  then  Fourier  transformed,  it  is  essentially  the 
same  one  found  for  the  pure— tone  response.  For  the  latter,  to 
obtain  the  displacement  of  the  basilar  membrane  it  is  necessary  to 
solve  a  nonlinear  eigenvalue  problem  for  the  local  inviscid 
wavenumber  k0  (Holmes  and  Cole,  1984;  Holmes,  1986).  This  problem 
is  given  by  the  following  partial  differential  equation  for  the 
(nondimensional ized)  displacement  of  the  basilar  membrane  in  the 
transverse  cross— sect i on 


(D2dy  ~  u2u)n0  =  -  2u2p0(x,y,O)  ,  (Al) 

where  D2=h3(x>  and  u=h(x).  The  function  p0(x,y,z)  comes  from 
the  fluid  pressure  and  is  given  as 

a> 

p0  <x ,y,z)  =  y  pm  coshxm(H-z )cosrmy  ,  (A2) 
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Inserting  (A2)  into  (Al)  the  problem  that  remains  to  be  solved  ii 

»  G 
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To  solve  (A4)  we  expand  u0 


in  beam  modes  as  follows 


n0  (x,y)  =  £ _ ,bn  (x  )  un  (x  ,  y )  .  (A5> 

n=  1 

The  functions  un(x)  are  the  natural  modes  of  the  simply  supported 
elastic  beam  in  the  transverse  cross-section  and  are  given  as 


un  <x  ty> 


cos(rny) 
y  G<x)  ■ 


<A6a) 


where 

(2n— 1 ) w 
rn  2G(x) 


Substituting  (A5)  into  (A4) ,  and  using  the  orthogonality 
modes,  the  following  system  of  equations  is  obtained 


( A6b ) 

of  the 
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In  addition  to  the  bj's  the  above  system  of  equations  also 
determines  the  local  inviscid  wavenumber.  In  fact,  the  dispersion 
relation  is  the  character i st i c  polynomial  for  (A7) . 

Up  to  this  point  all  of  the  beam  modes  have  been  included  in 


the  representat i on  of  the  solution.  However,  because  of  the 


symmetry  of  the  problem  (in  y)  the  first  beam  mode  is  all  that  is 
necessary  to  obtain  an  accurate  approximation  to  the  solution.  In 
this  case.  Me  have  from  <A5>  that 


H  o  (x  ,y )  = 


L  2G(x)J 


G  (x ) 


Also,  the  system  of  equations  (A7)  reduces  to  finding  the  solution 
to  the  following  dispersion  relation  for  the  local  inviscid 
wavenumber  k  0 


D2ov  “  “2u  =  )  , 


(A10) 


where 
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Once  k0  is  determined  from  (A10)  the  amplitude  Ar(x,u)  and  the  phase 
8r(x,u)  of  the  right-travel ing  wave  are  given  by 
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The  -function  kt(x)  in  (All),  (A12)  represents  the  viscous 
contribution  to  the  wave  number  and  is  given  as 
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The  parameters  are 

B  6  **B 
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and  H(x)=H*/B  is  the  nondimensional  "radius"  of  the  cochlear  fluid 
chamber  (see  Fig.  1).  The  constant  A*  in  (A14)  is  the 
nondimensional  area  of  the  stapes  footplate  and  it  is  related  to  the 
actual  area  Aj  as  follows 

Aw  =  B  Aw  . 

Finally,  in  (A13)  the  function  pm,  which  is  defined  in  (A3),  is 

—  —  ~acml<:m  i 
Pm  XmsinhxmH 

The  amplitude  A(x,u)  in  (25)  includes  Ar  as  well  as  the 
contri buti ons  of  the  waves  reflected  at  x=l.  To  account  for  these 


let  A; (x,u)  and  (x ,u)  be  such  that 


where 


x 

k  =  f-  J  [k0  +  S0  (1  -i)kt]  ds 

o 

and 

k,  =  k  I 

x  lx  =  i 

With  this 

A(x,u)  =  Ar  (x,u)  Ax<x,u)  , 

and 

6<x,u)  =  Or(x,u)  -  Ojt  (x,u) 


( A1B) 


( A19) 


